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Introduction 

Solid  Rocket  Motor  (SRM)  mission  design  requires  the  performance  prediction  of  the 
motor  until  burnout. 


The  experimentation  cost  is  very  high  due  to  manufacturing  of  SRM  and  due  to  the 
instrumentations  required.  Therefore  before  any  product  design  is  completed  for 
experimentation,  it  should  be  investigated  in  details  by  internal  aerodynamic  solutions 
for  minimum  design  correction. 


Time 
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Introduction 


The  internal  aerodynamic  solution  for  SRM  is  complicated  with  nearly  every  aspect  of 
Computational  Fluid  Dynamics’  (CFD)  topics.  The  full  SRM  solution  requires: 


1 .  High  Speed  Turbulence  Modelling 

2.  Combustion 

3.  Two-phase  Flows 

4.  Acoustics 


In  addition  to  complexity  of  physics, 
the  control  volume  in  which  the  flow 
solutions  will  be  carried  out  is  not 
constant. 


As  the  combustion  continues  the  grain  surface  moves  towards  the  propellant,  thus  increasing 
the  flow  volume.  The  CFD  solutions  require  unsteady  flow  solutions. 
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Introduction 


In  this  project  the  regression  of  combusting  surface  of  solid  propellant  in  rocket  motors 
is  investigated  by  using  a  fully  3 -dimensional  method  which  is  called  the  “Fast 
Marching  Method  (FMM)”. 

The  FMM  is  capable  of  following  the  surface  regression  under  local  rate  variations  by 
itself.  However  the  aim  is  to  use  a  varying  bum  rate  under  the  effects  of  pressure  and 
velocity.  Different  flow  solvers  are  coupled  with  the  FMM. 

At  each  time  step  the  FMM  will  generate  mesh  files  for  the  solver  and  obtain  fluid 
properties  and  therefore  the  bum  rate. 
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Surface  Regression 


The  fluid-propellant  interface  moves  in  the 
normal  direction  to  the  surface.  The  speed  of 
the  motion  is  dependent  upon  the  rate  of  bum 
of  the  propellant  which  is  dependent  upon  the 
properties  of  the  internal  flow. 


The  interface  can  be  located  either  by 
tracking  the  interface  (Lagrangian)  or  by 
capturing  (Eulerian)  on  stationary  mesh. 
The  FMM  method  is  of  Eulerian  type. 
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Surface  Regression 


Fast  Marching  Methods  (FMM)  are  numerical  algorithms  for  solving  the  nonlinear  Eikonal 
Equation  on  a  Cartesian  mesh. 


||Vw(jc)||  =  F(x)  inQ,  F(x)>0 
u  -  g(x)  on  r 


llV7,(x)ll 


in  Q,  rb  (x)  >  0 


T  -  0  on  r 


General  Eikonal  Equation 


Our  Form  of  Equation 


where  rb  is  the  bum  rate  and  T  is  the  time  of  arrival  of  the  burn  surface. 
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Surface  Regression 
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Interface  Capturing 


1 .  T0  is  not  inside  the  element 

2.  Td>T0>Tc>Tb>Ta  - 

3.  Td>Tc>T0>Tb>Ta  - 

4.  Td>Tc>Tb>T0>Ta  - 


The  interface  is  a  triangular  element. 
The  interface  is  a  quadrilateral  element. 
The  interface  is  a  triangular  element. 


D  D 
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Interface  Capturing 

There  are  five  choices  at  this  instant  how  this  captured  interface  will  be  used. 

1 .  Find  the  interface  area  vs  web  and  supply  this  data  to  0-dimensional  solver. 

2.  Find  the  interface  perimeter  at  stations  along  the  motor  vs  web  and  supply  this 
data  to  1 -dimensional  solver. 

3.  Capture  the  3 -dimensional  mesh  and  leave  the  interface  elements  as  they  are  and 
supply  this  mesh  to  3 -dimensional  hybrid  solver. 

4.  Capture  the  3 -dimensional  mesh  and  at  the  interface  insert  new  tetrahedron 
elements  and  supply  this  mesh  to  3 -dimensional  tetrahedron  solver. 

5.  Capture  the  3 -dimensional  mesh  and  at  the  interface  move  inner  nodes  to  the 
boundary  and  supply  this  mesh  to  3 -dimensional  tetrahedron  solver. 
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Flow  Solver 


3-D  FLOW  SOLVER 

•Euler,  Internal,  Finite  Volume  Solver 

•Hybrid  Mesh 

•Explicit,  RK4 

•2nd  order  accurate  in  space. 

•Calculation  of  Flux:  Upwind,  SW,  WL,  Roe 

•Boundary  Conditions:  Mass  flux  inlet,  pressure  exit,  wall, 
symmetry 
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Flow  Solver 
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Flow  Solver 
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Flow  Solver-3 
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Test  Case  -  1 
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Hybrid  Mesh 


_ 

Test  Case  -  1 


Sliced  elements  are  converted 
into  Tetrahedrons 


Moving  inner  nodes  to  interface 
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Test  Case  -  1 
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Pressure 
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Test  Case 


Time 
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The  results  are  obtained  by  constant 
propellant  bum  rate  assuming 
inviscid  ideal  gas  flow. 


